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ABSTRACT 

Cosmic ray acceleration through first-order Fermi acceleration in a collisionless plasma 
relies on efficient scattering off magnetic field fluctuations. Scattering is most efficient 
for magnetic field fluctuations with wavelengths on the order of the gyroradius of 
the particles. In order to determine the highest energy to which cosmic rays can be 
accelerated, it is important to understand the growth of the magnetic field on these 
large scales. We derive the growth rate of the long- wavelength fluctuations in the linear 
regime, using the kinetic equation coupled to Maxwell's equations for the background 
plasma. The instability, driven by the cosmic ray current, acts on large scales due to 
the stress tensor and efficient scattering on small scales, and operates for both left- and 
right circular polarisations. This long-wavelength instability is potentially important 
in determining the acceleration efficiency and maximum energy of cosmic rays around 
shock waves such as in supernova remnants. 
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1 INTRODUCTION 

The acceleration of Galactic cosmic rays up to energies 
of 10^* eV is wid ely believed to result from first order 
Fermi ac celeration jAxford et al.lll977l : lKrvmskiilll977l : iBelJ 
Il978al lbl: iBlandford fc Ostrikeilll978t ) in supernova remnant 
blast waves. This mechanism relies on a converging flow, 
where cosmic rays gain energy every time they cross the 
discontinuity, and is also referred to as diffusive shock ac- 
celeration (DSA). On either side of the shock the cosmic 
rays isotropise by scattering off magnetic fluctuations, while 
the plasma is highly collisionless. The length scale of the 
magnetic fluctuations and the amplitude determine the scat- 
tering efficiency. The more quickly the particles isotropise, 
the quicker they may diffuse back accross the shock and the 
shorter the acceleration time. 

There exists a variety of cosmic ray streaming instabil- 
ities that amplify the local magnetic field and increase the 
scat tering efficiency of the plasma. For a recent review, see 
e.g. iBvkov et al.l (|2011f ). On small scales, i.e. smaller than 
the gyroradius, the streaming cosmic rays induce a return 
current that creates a force perpendicular to the background 
magnetic field. This results in a non-resonant instability that 
can amplify the magneti c field to m any times its initial level 
in the nonlinear regime (|Bel|[200J). Since scattering is most 
efficient when the larmor radius is of the order of the wave- 
length of the magnetic fluctuations, it is therefore expected 



E-mail: k.schurel@physics.ox.ac.uk 



that this instability is effective in confining the low-energy 
cosmic rays. 

Resonant instabilities arise in the regime where the 
parallel component of the wavelength corre sponds t o the 
gyroradius of the streaming cosmic rays (ILerchd 119671 : 
iKulsrud fc Pearcdll969l : IWentzellll974l : ISkillinBlll975f r When 
the perturbed field becomes of the same order as the back- 
ground magnetic field {SB/B > 1) the resonance is lost. The 
saturation level of the waves therefore is rather low, and this 
instability is not effective in amplifying the magnetic field 
to many times its zeroth order level. 

Whether the highest-energy cosmic rays are efficiently 
confined in the shock region depends on whether instabilities 
develop on scales larger than the gyroradius of the generat- 
ing cosmic rays. In this paper we show that an instability 
acts in this regime that arises from the cosmic ray current 
and is mediated by the stress tensor. The coupling of the cos- 
mic rays to the plasma on the small scales, due to effe ctive 
scatt ering as a result of the non-resonant instability (|BelJ 
[2OOJ), further aids in making both left and right-handed 
polarisations unstable. 

We use the Boltzmann transport equation together with 
Maxwell's equations for the background plasma, to derive 
the growth rate of magnetic fluctuations on the large scale. 
Rather than limiting ourselves to the diffusion approxima- 
tion, we include higher order anisotropy and find that the 
stress tensor takes an active part in mediating the insta- 
bility. The small-scale scattering couples the two polarised 
modes and integrates the short-scale with the large-scale in- 
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stability. The growth rate is found to steeply depend on the 
wavenumber, resulting in an rapid decline with wavelength. 

We limit ourselves to a mono-energetic beam of parti- 
cles. This will allow us to get a cleaner picture of how exactly 
the instability works. The distribution function generally is 
expected to follow a power law beyond the thermal peak, and 
has a cut-off at the high-energy end. If diffusion is a function 
of momentum, which is generally assumed, the high-energy 
cosmic rays diffuse farther upstream of the shock than the 
lower energy ones. The low-energy cut-off therefore increases 
with distance upstream of the shock. Since the number of 
particles steeply declines with energy, it can be argued that 
the instability is dominated by the cosmic rays with the low- 
est energy in the local plasma, being a function of distance 
ahead of the shock. 

In Section [2] we show how we derive the dispersion rela- 
tion for the instability that results from including the stress 
tensor and scattering on small scales. In Section [3] we evalu- 
ate the result and look at the growth rate of the instability 
for a range of collision efficiencies and wave numbers. We 
show the differences between inclusion or neglect of both the 
stress tensor and the small-scale diffusivity in Section 2] We 
discuss consequences for confinement of cosmic rays in Sec- 
tion [S] In Section [6] we show that in most regimes it suffices 
to limit ourselves to the second order anisotropy (i.e. the 
stress tensor) and illustrate the difference when higher or- 
der anisotropy of the distribution function is taken into ac- 
count. Finally, we discuss its implications and conclude in 
Section [3 



2 METHOD 

The distribution of relativistic particles can be described by 
a modified Boltzmann equation, also known as the Vlasov- 
Fokker-Planck equation, which reads (in cgs): 

9tf + V ■ Vf + -^ ^E -h ^^^-^\ ■ V„f = V. ■ (D • V„f), (I) 
m \ c y 

with f the particle distribution in phase space, v the particle 

2 

velocity, and D = ^^ (I — nn) the diffusion tensor with v 
the collision frequency. Here I is the identity matrix, and n 
the unit vector in the direction of the corresponding tensor 
component. The effective collision frequency is a result of 
small angle scattering by fluctuations in the magnetic field. 
The distribution function can be decomposed into an 
isotropic part and arbitrarily many moments, which sample 
the anisotropy to increasing order, as follows: 



= Jo H 1 o 1 5 T 



(2) 



vi + vg + vi 



In spherical coordinates the directional components are 
given by: 



Vx = V cos 6 

Vy — V sin 9 sin <j 

Vz — V sin 9 cos < 



(3) 



By integrating the distribution function over angle — 
[0,n],(l) = [0, 27r] after multiplying with sin6', the scalar 



transport equation can be derived (where we use v = c for 
cosmic rays): 

c. 



4n [dtfo + ^V ■ fi] 



0. 



(4) 



The first moment of the distribution function can be found 
by mulitplying it with cosS, sinSsin^, sinScosi;/) to get the 
moment in x, y, z respectively. The equation for momentum 
transport is thus found to be: 



47r 



atfi + cVfo + — (B X fi) + ufi + |cV ■ {fa} = 
mc 5 



(5) 



We restrict ourselves to a parallel shock and choose the 
direction of shock propagation to be in x, meaning that the 
background magnetic field only has a Bj, component (Bq — 
i3o|x). We assume the electric field can be neglected. 

The unperturbed cosmic ray current follows directly 
from the momentum transport equation when looking at 
the background quantities only, where we use that jcr ~ 
ricrqus = ncrqfx{o)c/S. The subscript or 1 between brack- 
ets is used to indicate zeroth order and perturbed variables 
in case ambiguity can arise between orders in the expan- 
sion of the distribution function (no brackets) and orders in 
the perturbation (with brackets). A single index {x,y,z) is 
used for fi when a specific directional component is meant. 
Equation [5] gives for the zeroth order current in x: 



fx{0) 



vLx 



fo, 



(6) 



with Lx the scale height of the cosmic ray density gradient 
that is determined by the diffusion coefficient over the shock 
velocity D/us- Using that v — c^/(3D) the zeroth order 
current is given by fx{0) = (3tts/c)/o. 

By eliminating fi between the first two moment equa- 
tions, while neglecting higher order terms, the diffusion ap- 
proximation is derived. Often the second- and higher order 
anisotropic parts are neglected, since each subsequent or- 
der in anisotropy is another factor Us/c smaller, where Us 
is the shock velocity. In this section however, we relax that 
assumption and include the next order in anisotropy (fa). 
To higher order the transport equations can be derived by 
taking subsequent moments of the distribution function. 

In systems with strong shock waves, such as for exam- 
ple SNR blast waves, it has been shown that the first or- 
der anisotropy, being effectively the cosmic ray curre nt, is 
responsible for driving an instability on small scales (jBelJ 
;2004)- The cosmic ray current triggers a return current in 
the plasma. The jcr x Bi force acts to focus the cosmic rays, 
where subscript I indicates perturbed variables. The result- 
ing reaction force pushes out the plasma and magnetic field 
for right-hand polarised waves, resulting in a plasma insta- 
bility that can act to amplify the magnetic field to many 
times its background level. Here we evaluate how the cos- 
mic ray streaming instability develops when collisions and 
fa are explicitly taken into account. 

The first order perturbation of Equation [5] yields com- 
ponents for fy and f^ as follows: 

dtfy -h ^ (Bo X f, + B, X 4(0)) + viy + '^cdxfxy = (7) 

ftf^ + —{Bo xfy + By X 4(0)) + vfz + -cdxfxz = 0. 

Together with the MHD equations for the thermal back- 
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ground plasma this forms a closed system: 
dtB = V X (u X B) 

Fr. 



pd,u = Ji^ii^ _ vp - V ■ n 



(8) 
(9) 



In a thermal plasma, the stress tensor H is zero. We assume 
a cold plasma where the thermal pressure can be neglected, 
meaning that also VP = 0. The last term is a resistive 
force resulting from the collisions with the cosmic rays, and 
is equal but of opposite sign to the frictional force in the 
momentum equation for cosmic rays. 

The current in the plasma, minus the V x P component, 
needs to balance the cosmic ray current minus the thermal 
component, such that: 
c 

ith - -;— V X B = -jcr 

Writing the momentum equation for the plasma as its equiv- 
alent in terms of the cosmic rays gives: 



pdtu = Fr 



c 

jcr X B 



(10) 



+ 1— (V X B) X B + iypcr{Ucr - u) . 

c 47r 

where p^r is the mass density of the relativistic particles, and 
where Ucr — u is the drift speed of the cosmic rays relative 
to the plasma. To zeroth order, where uo = this equals 
the shock velocity Us. We ignore the second term on the 
right hand side in the second line further on. This Alfvenic 
term only becomes important on very short scales when the 
tension in the magnetic field lines damps the waves. The 
frictional force Fr has been set to be equal and opposite 
to the force induced by cosmic ray collisions on small scales, 
being: vpcrUcr, where we neglect the contribution of ui. The 
drift velocity is related to fi by Ucr/o = fic/3, such that we 
can write the momentum equation as: 

dtui = -^^(Bo X fi(i) + Bi X fi(o)) + ^t'nfi(i), (11) 

with n the number fraction of cosmic ray protons to back- 
ground protons, and where we use that Ws/o = /xc/3. We 
can eliminate the magnetic field from this equation using 
Eq.[71 



i9tui = -^^ ( dth + -cdxh{i 



(12) 



Feeding this into the induction equation, we can express the 
perturbed magnetic field in terms of fi and f2: 



dtBy 



dtB, 



nc 2 

-^Bodx{dtfy + -vdj^f^y) 

nc 2 

-—-Bodx{dtfz + -vdxfxz)- 
3 5 



(13) 



We assume plane wave solutions of the form ^x — 
^^^i(kx-ujt) ^Q derive the dispersion relation which gives for 
the perturbations in the magnetic field: 



_ nkcBo f r 2kc 
nkcBo f , 2kc 



(14) 



The linearised momentum transport equation then is: 
( ■ \f iBa , , lfx(o) r, , 2ikc 



, ■ \ t , qBo 
mc 



— — -Dy H —Jx 

mc 5 



0. 



The relevant fxy and fxz components are derived from 
the next moment of the transport equation, being the heat 
transport or pressure tensor equation: 
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9t{f2} + c( Vfl--V-fl{l2} 



(16) 



+ ^ (B X {f2}) + 3/.{f2} + JcV ■ {fa} 
mc 7 



0, 



where we eliminated t he dtfp compone nt using Eq. [l] We 
use the notation as in Ijohnstoni (|l96Cl ). where { } is used 
to distinguish tensors from vectors when both appear in an 
equation, {Ii} is the l-ih order unity tensor, and the sub- 
script in [...]; indicates that the equation is a summation of 
permutations for ijk for each element l\ ways, divided by l\. 
This gives a symmetric tensor so that e.g. fxy = fyx with 
components that look like: 



dtfxy + -{dxfy+dyfx) 

H {B^ fxx ~ B^ fxz + B^ fyz 



(17) 

B'fyy)+Siyfxy = 0. 



For now we ignore contributions to the distribution 
function that are third and higher order in anisotropy (we 
come back to the higher order equations in Sect.[6]l, and set 
those (which includes dtk) to zero. The non-zero contribu- 
tions from Eq. [16] then are limited to: 



ijxjy 



qBo 



fxz + Zvfxy = 



(18) 



'^ Pi f . iBq r I q f 

^Oxjz ~T~ Jxy 'r 'ji^Jxz 

2 m,c 



such that we can substitute the f^ components in the mo- 
mentum transport equation and in the perturbed magnetic 
field components using: 



Jxy — 
fxz ^= 



ikc{3ufy +UJgfz) 
ikc{?,Vfz -UJgfy) 

2((3i.)2-foj2) ' 



(19) 



where we wrote qBo/mc — u)g for the gyrofrequency. These 
components represent the transport of the perpendicular 
cosmic ray current. 

We now have all the components to write Eq. [15] as a 
function of just fy and fz as perturbed variables: 



i^ — ioj + 



k''c''?,v 



iQ? 



5((3i/)2-fa;|) u;^ \b{{?,vY + u^) 



k^C^UJa 



fy 



7 A 2 

k C OJq 



5((3i.)2+a;|; 



+ 



iVL^ 



k^c^Sf 



5(i3u)^+ujl 



(20) 



- ibJ -\- 



k^c^3u 



iQ.^ 



5((3i/)2-ha;|) 



k^C^UJa 



C UJ„ 



5((3u)^ + ul 



+ 



ifl^ 



k^c^3p 



5((3i.)2 -f a;|; 



where fl = \/fcjoPo/(pc). These can be readily combined to 
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find the dispersion relation: 

K C 



uj = ^Q. I iu) 



5(3i^ =F iLUg] 



(21) 



V — iid ^ iujg + 



;„2„2 

K C 



5(3i' =F ia;g 



The upper signs correspond to the left-handed polarisation 
(/z = ify) ^nd the lower signs to the right-hand polarisation. 



3 EVALUATION OF DISPERSION RELATION 

Because the growth time of the instability will always be 
larger than the cosmic ray crossing time (1/fcc), the second 
term on the right hand side of Ea. l21l alwavs dominates over 
the first one, and the dispersion relation can be simplified 
to: 

u2„2 



±0.' 



fc c 



5{3u =p iujg) 



V =F icjg + 



k'c' 



5(3!/ =F iojg 



,(22) 



In the short wavelength limit {kc S> uig,v), where the 
last term in the denominator becomes the dominant one, 
this method naturally retrieves the growth rate Q that cor- 
respo nds to the small-scale non-resonant instability (JBelll 
|2004| ). 

In the long-wavelength limit, i.e. for krg <^ 1, the dis- 
persion relation can be reduced to: 



±n^ 



e 



5(3i^^ =F i/iuOgV — a;| 



(23) 



From this it can be readily seen that below resonance the 
dependency of growth rate on k is very steep: 7 ex fc"^". This 
holds for both polarisations. For which one the growth rate 
is more rapid depends on the relative scattering efficiency. 

The depen dency that w e find here is steeper than the 
one derived in iBelJ (I2OOJ) in the long- wavelength limit, 
which was oc k. The reason for the faster growth found in 
that paper is the integration of the distribution function over 
momentum, with the upper limit of the momentum set to 
(X). This means that even for long wavelengths there was al- 
ways some contribution to the instability from cosmic rays 
that are in resonance with the magnetic instability. In real- 
ity the maximum energy is finite and therefore on the long 
scales there are no resonant particles to contribute to the 
instability. 

In Fig. [1] we plot the growth rate as a function of 
both the scattering efficiency and the wavenumber. The up- 
per panel shows the right-hand polarisation (lower signs in 
the equation) and the lower panel the left-hand one (upper 
signs). For wavenumbers below resonance the growth rate 
follows a k^'^ dependency for both polarisations. Beyond 
resonance (on smaller scales), the growth rate converges to 
the small-scale non-resonant growth rate Q, for the right- 
handed polarisation (7 oc k^ ), whereas a steep decline for 
the left-handed polarisation is observed. 

The relative strength of the growth rate for the two po- 
larisations on long scales depends on the relative scattering 
efficiency. Efficient scattering favours growth of the right- 
handed polarisation whereas low rates of v/ug favours the 
left-handed growth rate. This can be more easily seen in 
Fig. [2] where the growth rate is represented as a function of 
the scattering efficiency in the long wavelength limit. In this 




10 r 




Figure 1. Growth rate 7 as a function of scattering efficiency 
v/ujg and wavenumber k. The upper (lower) panel represents the 
right (left) handed polarisation. The growth rate is normalised in 
units of v1/{rgc). 



figure a value of kr-g = 0.001 was assumed, and the growth 
rate is normalised with respect to Q.. The scattering rate 
for which the right-handed polarisation is maximum, which 
is also the point where the left-handed polarisation drops 
below the right-handed one, is v/iOg = •\/l/3. This value 
remains the same for different k, as long as we evaluate it 
below the resonance frequency. 

We have plotted a range in scattering efficiencies be- 
tween u/ojg = [0.01, 10]. Low values mean that the cosmic 
rays complete more than a full gyration before scattering, 
and the particles are essentially magnetised. The more fre- 
quent collisions become, the more effective diffusion occurs 
across field lines. Very often Bohm diffusion is used, for 
which the scattering frequency is of the order of the gy- 
rofrequency. It is assumed that such a scattering efficiency is 
optimal in isotropising the cosmic rays and hence maximises 
the acceleration efficiency. Very high values of the scattering 
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Figure 2. Growth rate 7 as a function of scattering efficiency 
for the longest wavelength modes (normalized to Q). The solid 
(dashed) line represents the right (left) handed polarisation. The 
real part of the frequency tr equals the growth rate of the opposite 
polarisation in the regime where kvg ^ 1. 
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Figure 3. Growth rate (normalised to Q) as a function of 
wavenumber for various values of the scattering rate. The solid 
(dashed) lines represent the right (left) handed polarisation. 
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efficiency means that the particles essentially do not 'see' a 
coherent magnetic field. This situation might arise when tur- 
bulence is very strong and disordered on the relevant length 
scales. 

In Fig. [3] we plotted the growth rate, normalised to Q, 
as a function of wavenumber, for different values of ly/uJg. 
From this figure it can be seen that the growth rate of the 
right-hand mode at large k is equal to fl, and it depends 
on v/ojg how rapidly after resonance this convergence is 
reached. For low scattering efficiency the left-hand polar- 
isation dominates, and a strong peak arises around reso- 
nance that connects the left-hand mode to the right-hand 
one. Evaluation of f up to higher orders than f2 seems to 
reduce this peak and eventually completely smooth it out, 
as will be discussed in Sect. |6] 
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4 DISSECTING THE MEDIATORS OF THE 
INSTABILITY 

In order to separate the infiuences of {2 and the collision 
term on the instability, we look at how the instability evolves 
when small-scale collisions are absent or when the pressure 
tensor is omitted. 

The effect of neglecting collisions on short scales can be 
easily seen when setting u/cdg = in the dispersion relation. 
We plot this in Fig. |4l from which we can see that both 
the short-scale and the long-scale instabilities are still there. 
However, the coupling between the right-hand and left-hand 
mode is completely absent. In the long- wavelength limit, 
the right-hand mode is stable and only the left-hand mode 
arises, and it retains the same k^'^ dependency. The turning 
point is at resonance, after which the short-scale instability 
takes over, which is the same one as the 2004 Bell instability 
with a growth rate 7 oc fc^'^. 

Omission of the pressure tensor, while keeping the 
collision term, also results in an instability, but with a 



Figure 4. Growth rate as a function of wavenumber when col- 
lisions are absent {u = 0). The solid (dashed) line represents 
the right (left) handed polarisation. The right-hand polarisa- 
tion is stable on scales larger than the resonant wavelength 
{krg < 1) whereas the left-hand polarisation is stable on small 
scales {krg > 1). 



much smaller growth rate. This can be seen by evaluat- 
ing Eq. [21] to lowest order in k (where it only appears in 
Q — ^kjoBo/{pc)), while disregarding the fc^ terms that 



originate from the inclusion of f2 . 
dispersion relation reduces to: 



2 iSl Lu{u ± iujg) 
^ = T 5— — 9 ■ 



(24) 



When 1/ ^ u)g this would yield purely growing modes, 
where the growth rate (damping rate) of the right (left) hand 
polarisation equals: 



7; 



± — . 



(25) 
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Figure 5. Growth rate of the right-hand polarisation as a func- 
tion of scattering efficiency and wavenumber (normalised in units 
of v1/(rgC)). The black (upper) surface shows the growth rate as 
it was derived earlier including the stress tensor. The red (lower) 
surface shows the growth rate in case the stress tensor has been 
ignored. The dependency on wavenumber is less steep (7 oc k) 
but the growth rate remains much smaller nonetheless. Only the 
right-hand mode is unstable, so we have not plotted the left-hand 
mode. Neglect of the stress tensor results in missing vital parts 
of the current-driven instability. 



Only the right-hand mode is unstable. Although the depen- 
dency of the growth rate on k is less steep (7 oc fc) , its value 
is much lower compared to when f2 was included. This can 
also be seen in Fig. O in which the growth rate with and 
without including f2 is compared. 

In the limit where v <^u)g both modes are stable. In the 
absence of collisions (in addition to omission of the pressure 
tensor) the dispersion relation reduces to: 



n^ 



(26) 



and the system is stable. 

This illustrates that the inclusion of the stress tensor 
(or momentum transport) forms an integral part of the long 
wavelength instability. The strong dependency of u on fc in 
the dispersion relation (EQ. I2ip makes the k^ term dominant. 
This term arised from the addition of f2, the pressure tensor. 
Omission of this term will result in an unphysical result that 
neglects an important mediator of the instability. 



5 GROWTH OF THE INSTABILITY AROUND 
RESONANCE 

In the previous section we have seen that the diffusion term 
is needed to mediate the instability between the two polar- 
isations. On scales shorter than the resonance wavelength, 
the dominant mode is the right-hand polarisation. Growth 
in this regime is rapid, and it has been shown to amplify the 
magnetic field to many times the background field in the 
nonlinear stage. The filamentary structures that arise in the 



nonlinear regime grow as a result of the jxB force. Whether 
the largest growth around resonance arises due to the non- 
linear development of the rapid small-scale instability, or 
due to the linear growth directly resulting from the resonant 
streaming instability, remains uncertain. Small-scale struc- 
tures evolve to larger scale structures as was shown in iBelJ 
(2005). 

For scales larger than the resonance wavelength we have 
demonstrated that the left-hand polarisation dominates the 



instability as long as u/ujg 



Due to the longer co- 



herence time of the fi(i) x Bi force, the scattering of pos- 
itively charged cosmic rays is more efficient f or left-hand 
polarised magnetic fiuctuations (|Wentzellll974l ). Therefore, 
even though the growth rate is smaller, the resulting con- 
finement of cosmic rays may still be quite efficient for longer 
wavelengths, in effect for cosmic rays with higher energies 
than the dominant ones that drive the instability. 

The local peak in the growth rate at the resonant wave- 
length reinforces the imp ortance of the r esonant instability 
discus s ed for example by Lerche Jl967|) ; iKulsrud fc Feared 
(|l969l ): IWent^ (|l974l ): ISkillind (|l975l ). Simulations that can 
treat scales both shorter and longer than resonance will be 
needed to determine the relative contribution of the various 
instabilities. 



6 INCLUDING HIGHER ORDERS OF THE 
DISTRIBUTION FUNCTION 

In the derivation for the dispersion relation we evaluated 
the anisotropy of the distribution up to f2. In Section [J] we 
illustrated the difference in the dispersion relation when the 
pressure tensor anisotropy is neglected. Because of the major 
differences, it will be beneficial to determine what happens 
to the instability when we include anisotropy terms of higher 
order moments. 

In this section we evaluate the dispersion relation for 
the next order in f (fs) and show that no major differences 
arise in comparison to neglecting it. The higher order trans- 
port equations are derived as described in Sect. [21 The third 
moment yields the pressure transport or heat tensor equa- 
tion: 



Stt 
35 



dt{h} + c{V{h} - ^V ■ {f2}{l3}) 



(27) 



+ ^(B X {fa}) + &u{h} + jcV ■ {U] 
mc 9 



= 0. 



In the pressure tensor equation (Eq. I16[l we now have 
to take into account the contributions from fa and the time 
derivative in f2, being of the same order as fa. The contri- 
butions from the heat transport equation then are: 

{3U - iuj)fxy - OJgfxz + -^fxxy + — /y = (28) 

{2,V - iuj)fxz + l^gfxy H :=— fxxz + "^ /z = 0. 

As can be seen from evaluating fxy and fxz , where now 
we do include components of fa, we only need to determine 
fxxy and fx^z in order to get the dispersion relation. These 
can be derived to be: 



o 
Otjxxy 4" TT^COx Jxy ^g Jxxz ~r '^^ Jxxy — U 



(29) 



A long-wavelength instability involving the stress tensor 7 



dtfxxz + -rrCdxfxz + ^gfxxz + Gl^fxxz — 0, 
15 

where we neglect terms that are of higher order than fs, 
being its dt component and f^. Eqs. 1291 are used to solve for 
fxy and fxz, which can be substituted in the equation for fi 
to give the modified dispersion relation: 



(30) 



T!^ 



1.2 „2 
K C 



5 [3u ~ iu! ^ iujg + 



8fe2c2 



35(6i'^iujg) 




ioj^OJg + 



K c 



5 {3u — iu! ^ iujg + 



8fe2c2 



35(6i^^iajg ) 



Fig. |6] shows the growth rate for the dispersion relation 
up to fa as a function of both v/cOg and k. On the long scales 
there is no noticeable difference between this result and the 
one ignoring fa. The peak around resonance softens a bit and 
moves to slightly shorter wave numbers, closer to kvg = 1. 
We have evaluated the dominant terms for subsequent orders 
in f (up to fio) and found that this trend seems to continue 
when higher orders of f are included. If the trend continues 
as we expect, the transition around resonance smoothens 
out completely. 

A curious effect arises for the left-hand mode in the 
short- wavelength limit. Rather than steeply declining with 
wave number, it increases when the highest order of f that 
is included is odd (i.e. fa or fs). When the highest order 
of f that is included is an even number, the growth rate 
asymptotically decreases with wave number. From the dis- 
persion relation it can be seen how this effect arises: When 
the highest order is an even number, the orders of k in the 
nominator are exactly matched in the denominator and thus 
cancel, whereas this is not the case when the highest o rder 
is an odd number. In the analysis bv lAchterberg ([1983) the 
growth rate in this regime is shown to decline for high k, 
which leads us to believe that the real trend is captured 
when the highest order is an even order of f (f2, f4), which 
has a decreasing growth rate with wave number for the left- 
hand mode. 





Figure 6. Growth rate for the dispersion relation up to orders of 
fa as a function of scattering efficiency and wavenumber. The up- 
per (lower) panel represents the right (left) handed polarisation. 
The growth rate is normalised in units of v'^/{rgc). 



7 DISCUSSION AND CONCLUSIONS 

We have derived the dispersion relation for an instability 
that is driven by the zeroth order cosmic ray current. The 
long-wavelength component is mediated by the stress tensor, 
and coupling between the left- and right-hand modes occurs 
when scattering of cosmic rays off small scale turbulence is 
taken into account. 

The streaming of cosmic rays ahead of the shock induces 
a return current in the plasma that triggers instabilities from 
large scales down to the Alfvenic regime on very small scales 
at which the magnetic tension becomes efficient in damping 
the oscillations. We did not include the Alfven term (mag- 
netic tension), which only becomes of importance beyond 
krg > 10^ for our parameters. 

On small scales our analysis reproduces t he hydrody - 
namical non-resonant instability as described in lBelll l|2004 ). 
The growth rate at small scales is very high and the am- 
plified magnetic field can reach many times its zeroth order 



strength. In an environment such as that upstream of a SNR 
blast wave, the efficient generation of magnetic field ampli- 
fication on small scales is expected to increase the scatter- 
ing efficiency of the lower energy cosmic rays. We find that 
this in turn seems to couple the left- and right-hand modes, 
resulting in a small left-handedly polarised contribution to 
the short-scale instability, in addition to the dominant right- 
handed mode. 

On long scales, predominantly the left-hand mode is 
unstable. This effect only arises when anisotropy of order 
f2 (and higher) is taken into account. The long-wavelength 
mode behaviour already converges when orders of f2 are in- 
cluded, as we find that higher-order anisotropy terms do not 
alter the instability significantly. Inclusion of the collision 
term results in coupling between the two polarisations such 
that both modes are also unstable in the long-wavelength 
limit. The growth rate of the long- wavelength part of the 
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instability is a function of fc^'^, which results in a rapid de- 
cline of the growth rate with wavelength. In a recent paper, 
iBvkov et alj ([2011]) found a larger growth rate in this regime 
using a different method that includes the ponderomotive 
force resulting from the instability on short scales. 

Since we are assuming a mono-energetic beam of par- 
ticles not all the effects around the resonance points are 
captured in this approach. However, since the number of 
cosmic rays steeply declines with energy, it can be expected 
that the instability is dominated by the cosmic rays with the 
lowest energy in the local plasma, being a function of dis- 
tance ahead of the shock. Including higher orders of f shows 
converging behaviour in the resonance region, and the peaks 
that arise as a result of just including up to f2 smooth out. 

Alternatively, the magnetic fluctuations in the reso- 
nance regime could be dominated by filaments resulting from 
the short-scale nonresonant instability that have grown in 
structure to length scales that are in resonance or beyond. 
It will be important to evaluate the non-linear regime and 
saturation levels to determine the actual efficiency of this 
instability in confining high energy cosmic rays. 

What happens farther upstream, in the region where 
only the high-energy cosmic rays diffuse out to, while the 
lower-energy cosmic rays are confined closer to the shock 
region, is uncertain. If effective diffusivity only operates 
when a significant fraction of the low-energy cosmic rays is 
around, the two regimes may not couple, causing the long- 
wavelength instability to occur exclusively in the left-handed 
variety and vice versa for the short-scale instability. 

An interesting result is that the dominant mode on 
scales larger than the gyroradius is the left-hand polarised 
one. This polarisation is more efficient in scattering cosmic 
rays since the perturbed j x B force continues to act in the 
same direction over longer scales. The deflection therefore is 
larger, and scattering more efficient. Therefore, even though 
the growth rate of this mode is smaller, the long-wavelength 
instability may still be very effective in confining cosmic rays 
that have gyroradii larger than the driving cosmic rays. 
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